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A nearly-free-electron (NFE) model to describe STM spectroscopy of (111) metal surfaces with 
Kondo impurities is presented. Surface states are found to play an important role giving a larger 
contribution to the conductance in the case of Cu(lll) and Au(lll) than Ag(lll) surfaces. This 
difference arises from the farther extension of the Ag(lll) surface state into the substrate. The 
different line shapes observed when Co is adsorbed on different substrates can be explained from 
the position of the surface band onset relative to the Fermi energy. The lateral dependence of the 
line shape amplitude is found to be bulk-like for ^ 4 A and surface- like at larger distances, in 
agreement with experimental data. 
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I When a magnetic impurity is inside a metallic host, 
the Kondo effect can occur 1]. Signatures of the Kondo 

[ effect also appear in STM measurements of noble metal 
surfaces with adsorbed Kondo impurities through the ap- 

' pearance of characteristic zero bias line shapes 0, 0| of 
the Fano-type 4]. These line shapes are found to de- 
pend strongly on the metal surface on which the mag- 
netic atom is adsorbed. For instance, the line shape as- 
sociated with Co on Au(lll) is more asymmetric than on 
Cu(lll) whereas for Co on Ag(lll) a rather symmetric 
line shape is observed. As both bulk and surface states 
are present in (111) surfaces it is yet an open question 

I which of these play the most important role in the Kondo 
effect and what determines the line shapes observed for 

I the different substrates. 

I In spite of the large amount of experimental work ded- 
icated to the characterization of line shapes associated 

' with different noble metal surfaces a complete theoreti- 
cal model for the tip-surface-adsorbate interaction is yet 

] lacking. The interaction of bulk states in the metal with 
the tip and adsorbate has been done through a Jellium 
model for the surface with step ^,^] , image 7] and Jones, 
Jennings and Jepsen (JJJ) [8| potentials. Tight-binding 
approaches have also been introduced^- However, the 
contribution of surface states to the conductance has only 
been partially discussed 

Previously, we studied the tip-substrate-adsorbate in- 
teraction considering bulk states only In this Letter, 
we introduce a NFE description of the substrate which is 
the simplest way to describe bulk and surface states on 
equal footing. Surface states are found to give a substan- 
tial contribution to the STM conductance even when the 
tip is right above the Kondo impurity. This contribution 
varies from surface to surface and is larger for Cu(lll) 
and Au(lll) than for Ag(lll). This can be understood 



from the fact that the surface state in Ag(lll) extends 
far into the substrate. Finally we find that the position 
of the surface band onset relative to the Fermi energy, 
ep, determines the different line shapes observed in the 
different substrates. 

A NFE description of the surface is introduced. The 
ionic potential in the (111) direction is taken into ac- 
count as a perturbation to the Jellium potential inside 
the crystal. This potential opens up a gap and surface 
states split from the bulk band edges. Since the weight 
of the surface states is removed from the bulk states one 
might expect that they are not important for integrated 
properties. However, they can be relevant to Fermi en- 
ergy properties such as the STM conductance. Apart 
from including surface states we make similar assump- 
tions as in Ref. (i) neglect the direct coupling 
of the tip with the substrate d bands and with the 3d 
orbital of the adsorbate due to the localized nature of 
the d orbitals [ij] and the large tip-surface separation, 
(ii) the adsorbate is modelled by a single c?322_r2-orbital 
and the tip by a single s-orbital ^12. 13], and (iii) the 
momentum dependence (including orthogonalization ef- 
fects) of the hybridization matrix elements are explicitly 
taken into account 10, 14]. 

When Co is deposited on a noble metal surface, it cap- 
tures electronic charge so that there is effectively only 
one hole left |l5l ] . The unpaired spin residing in Co is re- 
sponsible for the Kondo effect experimentally observed. 
Hence, we introduce a generalized Anderson model to de- 
scribe the substrate, adsorbate and the tip. Neglecting 
the orbital degeneracy of the 3d orbital the model reads 

k,(T o" 

+ E ^k(4ck, + H.c.) + Ud\d^dld^ 

k,cr 



2 



k,cr 



+ H.c.) + Ht 



(1) 



Here, is the energy level of the adsorbate d orbital, 
ct creates an electron with spin a and momentum k. dl. 
and create an electron in the d3^2_^2 orbital of the ad- 
sorbate and the s-orbital of the tip, respectively, eg. and 
are the metallic energies and the hybridization ma- 
trix elements between the substrate and the adsorbate, 
respectively. U is the Coulomb repulsion of two electrons 
in the adsorbate. The last two terms in the model are the 
tip-substrate interaction which is governed by the matrix 
elements, Mg., and Ht that describes the tip which is sup- 
posed to have an unstructured density of states. 

The modification of the STM conductance through the 
surface due to the presence of the 3d impurity is given 

by 113 

5G{uj) = GoPtlui{{A{uj)+iB{uj))Gdd{uj){A*{uj)+iB*{u;))}, 

(2) 

with Go = ^ and pt the tip density of states. 
In the above equation, B{uj), reads 
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(3) 

where a sum over the two bulk bands appearing due to 
the ionic potential in the (lll)-direction is understood in 
the k sum. 

The first term in Eq. Q describes the interaction of 
the adsorbate and tip with bulk states and the second 
describes the interaction with the surface band. A{uj) 
is the Kramers-Kronig transformation of B{uj). Matrix 
elements, Mg. and V^, are evaluated with the orthogonal- 
ized wavefunctions, |k >. Details are given in Ref. [loj . 
For the systems of interest here, A{u!) and B{llj) are real. 
The crucial quantity is B{uj) as it embodies the complete 
information concerning the tip-substrate-adsorbate sys- 
tem. B(uj) depends on the tip position, R = (R||,Zt), 
the adsorbate position, (0, Zd), where Z is refered to the 
last plane of ions. The metal potential, Vm, comes in 
the matrix elements, and Afg.. We define throughout 
the paper _R|| = |R|||. The Green's function, Gddi^), de- 
scribes the electronic properties of the Sd adsorbate im- 
mersed in the metallic continuum including many-body 
effects such as the Kondo effect. The Kondo peak is sim- 
ulated through a Lorentzian of width Tk positioned at 
eK and is contained in Gdd{^) with Tk and e^c taken 
from experimental data. These parameters vary between 
ex ~ 3 - 6 meV and Tif ~ 50 - 90 K for the different 
(III) surfaces 0, [Hill- 
Metal wavefunctions are obtained by solving 
Schrodingers equation in the presence of the sur- 
face potential, which is given (in Rydberg energy units) 



TABLE I: Surface state data and parameters used in a NFE 
model of (111) noble metal surfaces. The depth of the crystal 
potential, Vb, the first Fourier component of the ion poten tial, 
Vg, and the metal work function, W , are taken from Ref. 
and are given in eV's. Vb is referred to the vacuum level while 
the surface state energy is referred to the Fermi energy, ep- 
The lattice parameter, a, and the decay of the surface state 
into the bulk. Ass, are given in A while the steepness of the 
potential denoted as A, is given in A~^. 
Surface a Vp Vg W A Zjm ess — (-f mss /me Xss 
Cu(lll) 3.64 14.7 2.55 4.95 2.21 1.29 -0.32 038 ToT 
Ag(lll) 4.08 15.1 2.1 4.74 2.23 1.28 -0.05 0.39 28.6 
Au(lll) 4.07 17 2.3 5.30 2.36 1.25 -0.49 0.26 13.5 



by: 



Vm{Z) = { 



2{Z-Z.. 

-Vo 



I Zm < Z < Zi, 



-Vo + 2Vgcos{GZ) 



Z < Z„ 



(4) 

where Z is the perpendicular direction to the surface (the 
(III) direction). Vq is the depth of the bulk potential and 
A controls the sharpness of the surface barrier potential. 
The parameters A = 2Vb/A — I and f3 = Vq/A are ob- 
tained by imposing the condition of continuous differen- 
tiability of the JJJ potential at Z = Zim (see Ref. for 
more details). Zm is obtained by matching the periodic 
part of the potential in the perpendicular direction to the 
surface with the JJJ potential. All distances are referred 
to the last plane of ions and energies are referred to the 
vacuum level unless otherwise stated. G is the reciprocal 
vector along the perpendicular direction to the surface 
(the (III) direction) and Vq is the first Fourier compo- 
nent of the crystal potential associated with G. For the 
(III) surfaces G = 27r-\/3/a, where a is the lattice pa- 
rameter. 

The periodic part of the potential inside the crystal 
couples the |k > and |k -I- G > states, leading to a gap of 
size: Eg = 2|Vg|. For perpendicular momenta: < fc^ < 
G/2 and energies outside the gap, bulk wavefunctions 
are obtained by matching the wavefunctions outside and 
inside the crystal. Wavefunctions are free electron-like 
far from the band edges but differ from this description 
close to them. 

Surface states are possible inside the gap with momen- 
tum kz = G/2 + i/Xss- These states decay inside the bulk 
with the decay length. Ass , retaining the periodicity par- 
allel to the surface. By matching the complex solution 
inside the bulk with the solution outside we solve for the 
energy at which the surface state exists jioj . 

Wavefunctions associated with these surface states 
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FIG. 1: Conductance line shapes obtained from the NFE 
model for Co on Cu(lll). In (a) we show conductance line 
shapes for different lateral positions of the tip associated with 
the function B{uj) displayed in (b). The inset shows the 
lateral variation of the normalized conductance amplitude, 
yl(7?ll)/Ao compared to experimental data. Distances are 
given in A. 

read 

{^e*''ii''iie(^-^'")/^-cos(fz-h 61(e)) , z < Z, 
e'''ii'-|li/;(e,z) , z> Z, 

(5) 

where, '0(e, z), is the solution to the Schrodinger equa- 
tion in the perpendicular direction for a given energy, e. 
a is the area of the square enclosing the surface wave- 
function and 9{e) is the phase factor coming from the so- 
lution of Schrodingers equation inside the bulk. Surface 
states with larger Xss penetrate more into the solid and 
correspond to surface state energies closer to the bulk 
band edge. The normalization constant, D ~ y^l/Xss 
for Xss ^ 1, going to zero as the surface state energy 
approaches the bulk continuum. This constant is im- 
portant to determine the relative contribution of surface 
states compared to bulk states in the conductance. 

The surface state energy is determined by keeping all 
parameters describing the surface potential fixed except 
for Zim which is varied to recover the surface state po- 



sition in agreement with photoemission data 'l8'|. The 
values of Zim obtained are shown in Table For in- 
stance, for Cu(lll) a surface state at 0.32 eV appears 
below the Fermi energy and the first image state at 0.83 
eV below the vacuum level for Zim = 1-29 A in good 
agreement with experimental values. Bulk-band effec- 
tive masses, m* /rrie = 0.74, 0.55, and 0.52 are taken in 
the (111) direction and 0.25, 0.36 and 0.18 in the parallel 
direction for Cu, Ag and Au, respectively. In our calcula- 
tions of B{lu) we neglect the image state band as it is far 
from the Fermi energy and disperses over a rather small 
energy range, leading to a small contribution to A{€f)- 

Parameters describing surface bands of Cu(lll), 
Au(lll) and Ag(lll) surfaces are summarized in Ta- 
ble HI where relevant parameters of the surface potential 
together with values of Xss and surface band effective 
masses, russ, are displayed. These masses are in good 
agreement with available experimental data from STM 
of clean surfaces 0,113,01 and photoemission data [l8j |. 

The conductance for the three noble metal surfaces 
with Co adsorbed on them is calculated. Results for B{uj) 
and conductance line shapes are displayed in Fig. ^ to 
Fig. 13 for Cu, Au and Ag, respectively. Most of the con- 
tribution to B{uj) close to the Fermi energy comes from 
the surface band and is larger for Cu(lll) and Au(lll) 
than for Ag(lll) as expected as in this case the surface 
state penetrates more into the crystal. At the same time 
B{uj) has a rapid drop with increasing to. This is because 
as the surface band disperses it gradually approaches the 
bulk continuum so that at a certain energy the surface 
band hits the bulk band and its amplitude drops to zero. 

Depending on the position of the surface band onset 
relative to the Fermi energy, the conductance line shape 
becomes more or less symmetric. If there is more weight 
below the Fermi energy, the conductance line shape is 
asymmetric with positive Fano parameter, q > (see 
Fig. 1 of Ref.[l3). From Table. Qwe see that \ess - epl 
is smallest for Ag(lll) and largest for Au(lll). Hence, 
the most asymmetric line shape corresponds to Au(lll) 
and the most symmetric one to Ag(lll). The different 
line shapes calculated for the different substrates are in 
agreement with experimental observations. 

The conductance amplitude, A(i?||), is more rapidly 
suppressed at short distances ^ 4 A than for > 4 
A. This is because both bulk and surface states contribute 
to B{eF). As the tip is displaced laterally the bulk con- 
tribution to B{ll!) close to drops quickly until it even- 
tually vanishes. Farther apart from the impurity, the sur- 
face band dominates and the lateral dependence becomes 
surface-like. Had we considered bulk states only in our 
model, >l(i?|| — 6)/Ao ~ 3%, which is an order of magni- 
tude smaller than the one found experimentally^^- This 
is because A{R\\) ^ 1/^^ for bulk states. In contrast, the 
amplitude dependence on is in good agreement with 
our calculations when both surface and bulk states are 
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FIG. 2: Conductance line shapes calculated from the NFE 
model for Co on Au(lll). Line shapes become more asym- 
metric with more positive Fano parameter as e^ — ess increases 
(compare the plotted line shape for Au(lll) with the ones 
shown in Figs. Qand|HJl. 



included (see inset of Fig. . 

In conclusion, surface bands give an important contri- 
bution to the conductance in STM measurements of (111) 
noble metal surfaces with Kondo adsorbates. The posi- 
tion of the surface band onset relative to the Fermi en- 
ergy determines the line shapes in the different (111) sub- 
strates. The lateral variation of the amplitude is found to 
be bulk-like close to the adsorbate becoming surface-like 
for > 4 A in agreement with experimental observa- 
tions. 
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